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Abstract 

Random walks and polygons are used to model polymers. In this pa- 
per we consider the extension of writhe, self-linking number and linking 
number to open chains. We then study the average writhe, self-linking and 
linking number of random walks and polygons over the space of configura- 
tions as a function of their length. We show that the mean squared linking 
number, the mean squared writhe and the mean squared self-linking num- 
ber of oriented uniform random walks or polygons of length n, in a convex 
confined space, are of the form 0(n 2 ). Moreover, for a fixed simple closed 
curve in a convex confined space, we prove that the mean absolute value 
of the linking number between this curve and a uniform random walk or 
polygon of n edges is of the form 0(^/n). Our numerical studies confirm 
those results. They also indicate that the mean absolute linking number 
between any two oriented uniform random walks or polygons, of n edges 
each, is of the form 0(n). Equilateral random walks and polygons are 
used to model polymers in ^-conditions. We use numerical simulations to 
investigate how the self-linking and linking number of equilateral random 
walks scale with their length. 

1 Introduction 

A polymer melt may consist of ring polymers (closed chains), linear polymers 
(open chains), or a mixed collection of ring and linear polymers. Polymer chains 
are long flexible molecules that impose spatial constraints on each other because 
they cannot intersect (de Gennes 1979, Rubinstein and Colby 2003). These spa- 
tial constraints, called entanglements, affect the conformation and motion of the 
chains in a polymer melt and have been studied using different models of en- 
tanglement effects in polymers (Orlandini and Whittington 2007, Tzoumanekas 
and Theodorou 2006). However, a clear expression of exactly what one means 
by entanglement and how one can quantify the extent of its presence remains 
an elusive goal. 
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As open polymer chains can always be moved apart, at least at long time 
scales, the definition of knotting or linking must refer to spatially fixed config- 
urations, and, as a consequence, does not give rise to a topologically invariant 
concept. In contrast, for closed polymer chains, the concepts of knotting and 
linking are unchanged under continuous deformations that do not allow breakage 
of the chains or passage of one portion of a chain through another. To charac- 
terize the knotting of an open chain one can use the DMS method (Millett et al 
2004, Millett and Sheldon 2005) to determine the spectrum of knotting arising 
from the distribution of knot types created by closure of the open polygon to 
the "sphere at infinity" . One is then able to employ the methods of traditional 
knot or link theory (Kauffman 2001) to study the topology of individual con- 
stituents or of the entire collection. The knot types of individual closures or 
of ring polymer chains can be analyzed, for example, using the HOMFLYPT 
polynomial, which can distinguish the different knot and link types with high 
precision (Freyd et al 1985, Przytycki and Traczyk 1987, Ewing and Millett 
1997). For pairs of oriented ring polymers, for "frozen" open polymer chains, 
or for a mixed frozen pair, we will see that the Gauss linking integral can be 
applied to calculate a linking number, which is a classical topological integer 
invariant in the case of pairs of ring polymers. For open or mixed pairs, the cal- 
culated quantity is a real number that is characteristic of the conformation but 
changes continuously under continuous deformations of the constituent chains. 
As polymers in solution are flexible moving molecular conformations, we will be 
interested in characterizations of their physical quantities at fixed moments in 
time as well as their time averaged properties. For the later, we will study the 
averages of such properties over the entire space of possible configurations. 

By Diao et al (1993), Diao (1995), Pippcnger (1989), Sumners and Whit- 
tington (1988) we know that the probability that a polygon or open chain with 
n edges in the cubic lattice and 3-space is unknotted goes to zero as n goes to 
infinity. This result confirms the Frisch-Wasserman-Delbruck conjecture that 
long ring polymers in dilute solution will be knotted with high probability. A 
stronger theorem is that the probability that any specific knot type appears as a 
summand of a random walk or polygon goes to 1 as n goes to infinity. However 
the probability of forming a knot with a given type goes to zero as n goes to 
infinity (Whittington 1992). 

The knot probability of polymer molecules also depends on the extent to 
which the molecule is geometrically confined. This has been studied by Ar- 
suaga et al (2007), Tesi et al (1994). For instance, DNA molecules confined to 
viral capsids have a much higher probability of being knotted. Moreover, the 
distribution of knot types is different from the distribution of DNA in solution 
(Arsuaga et al 2002, Weber et al 2006). 

The main subject of our study concerns uniform random chains confined 
in a symmetric convex region in M 3 . Each edge of the random open chain or 
polygon is defined by a pair of random points in the convex region with respect 
to the uniform distribution. Throughout this paper we will refer to open uniform 
random chains and to closed uniform random chains as uniform random walks 
and uniform random polygons respectively. We will consider uniform random 
walks or polygons confined in the unit cube. The uniform random polygon 
(URP) model (Millett 2000) is used to investigate the complexity of knots formed 
by polymer chains in confined volumes. 
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Figure 1: (a) +1 crossing and (b) —1 crossing 



This paper is organized as follows. In section 2 we study the scaling of the 
mean squared writhe, the mean squared linking number and the mean squared 
self-linking number of oriented uniform random walks and polygons in confined 
space with respect to their length. Next, we study the scaling of the mean 
absolute linking number of an oriented uniform random walk or polygon and a 
simple closed curve, both contained in a unit cube, with respect to the length of 
the random walk or the polygon respectively. In section 3 we present the results 
of our numerical simulations which confirm the analytical results presented in 
section 3. Although theoretical results about the absolute linking number be- 
tween two uniform random walks or polygons in confined space appear difficult 
to acquire, we are able to provide numerical results in section 3. Also, it is more 
difficult to provide analytical results for equilateral random walks or polygons 
than for uniform random walks or polygons. An ideal chain is an equilateral 
random walk composed of freely jointed segments of equal length (equilateral) 
in which the individual segments have no thickness (Diao et al 2003, Diao et 
al 2005, Dobay et al 2003) . This situation is never completely realized for 
real chains, but there are several types of polymeric systems with nearly ideal 
chains. At a special intermediate temperature, called the ^-temperature, chains 
are in nearly ideal conformations, because the attractive and repulsive parts of 
monomer-monomer interactions cancel each other. In this direction, we give 
in section 3 numerical estimations of the scaling of the mean absolute linking 
number and the mean absolute self-linking number of equilateral random walks. 

2 Measures of entanglement of open curves 

In a generic orthogonal projection of two oriented polygonal chains, each crossing 
is of one of the types shown in Figure [TJ By convention, we assign +1 to a 
crossing of the first type and —1 to a crossing of the second type. 

For a generic projection of two oriented curves Zi,?2 to a plane defined by 
a vector ( £ S 2 the linking number of a diagram, denoted Z&^Zi,^), is equal 
to one half the algebraic sum of crossings between the projected curves. The 
linking number of two oriented curves then is equal to the average linking 
number of a diagram over all possible projection directions, i.e. L(^i,^) = 
l/47r(/^ gS2 lk^(li,l2)dS). This can be expressed by the Gauss linking integral 
for two oriented curves. 

Definition 2.1. The Gauss linking number of two oriented curves l\ and I2, 
whose arc- length parametrization is 71 (t) , 72 (s) respectively, is defined as a dou- 
ble integral over li and I2 (Gauss 1877): 



L{h,h) = ]- [ I 



(71ft), 72(a), 7iW ~ 72(a)) 

\ll{t)-l2{sf 



dtds 



(1) 
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where (7i(i),72(s),7i(i) — 72 (s)) is the triple product of 71 (t), 72(3) and 

71 W -72O). 



Similarly, for the generic orthogonal projection of one oriented curve I to 
a plane defined by a vector ( e S 2 we define the writhe of a diagram, de- 
noted Wr^(l) to be equal to the algebraic sum of crossings of the projec- 
tion of the curve with itself. Then the writhe of a curve is defined as the 
average writhe of a diagram of the curve over all possible projections, i.e. 
Wr{l) = l/4:TT(J^ eS2 Wr^(l)dS). This can be expressed as the Gauss linking 
integral over one curve. 

Definition 2.2. The writhe of an oriented curve I, whose arc- length parametriza- 
tion is 7(i), is defined by the Gauss linking integral over a curve 



2tt y [0 ,i]. y [0 ,i]. i7(*)-7(*)i 3 

where [0,1]* x [0, 1]* = {(x, y) e [0,1] x [0, l]|z ^ y}. 



We define the average crossing number of a curve I, whose parametrization 
is 7(7), to be the average sum of crossings in a generic orthogonal projection 
over all possible projection directions. It is expressed by a double integral over 
I: 

ACN=±[ { 1(7(*).7W,7(0 -7(^)1^ (3) 
^ J [0,1] J [as] |7(*)-7(s)l 
where (^{t), j(s), -f(t) — 7(5)) is the triple product of ^(t), j(s) and j(t) — j(s). 



We observe that the geometrical meaning of the linking number and the 
writhe is the same for open or closed curves. The linking number of two curves 
is the average over all possible projection directions of half the algebraic sum of 
crossings between the two components in the projection of the curves. Similarly 
the writhe of a curve is the average over all projection directions of the algebraic 
sum of crossings in the projection of the curve. 

In the case of oriented closed curves the linking number is an integer topolog- 
ical invariant, i.e. it is invariant under isotopic moves of the curves. But in the 
case of oriented open curves the linking number is not topological invariant and 
it is not an integer. If the open curves are allowed to move continuously without 
intersecting each other or themselves in space, all the above measures are con- 
tinuous functions in the space of configurations. Furthermore, as the endpoints 
of the curves move towards coincidence, the linking number, self-linking number 
or writhe tend to the values of those measures for the resulting closed knots or 
links. 

Polymer chains are often modeled using open or closed polygonal curves. 
There exist several models of random walks or polygons that can be used and 
that are representative of the properties of different polymer melts. Since a 
polymer melt can take different formations in space through time, we are inter- 
ested in a measure of complexity of polymer chains that will be independent of 
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a specific configuration, and characteristic of the configuration space and how 
it depends on the length of the chains. In this paper, we will focus our study on 
uniform random walks in a confined space as this provides a simplified model 
for our theoretical study and will have a behaviour similar to other polymer 
models. 

3 Uniform random walks and polygons in a con- 
fined space 

The uniform random walks and polygons are modeled after the URP model, 
introduced by Millett (2000). In this model there are no fixed bond lengths and 
each coordinate of a vertex of the uniform random polygon contained in C 3 , 
where C = [0, 1] is drawn from a uniform distribution over [0,1]. 
The following theorem has been proved by Arsuaga et al (2007). 

Theorem 3.1. The mean squared linking number between two oriented uniform 
random polygons X and Y ofn edges, each contained in C 3 , is \n 2 q where q > 0. 
A similar result holds if C 3 is replaced by a symmetric convex set in M. 3 . 

The above result is independent of the orientation of the two uniform random 
polygons. Due to the weight squaring gives to larger linking number, we propose 
that the mean of the absolute value of the linking number between oriented 
uniform random walks or polygons would be a more informative measure of the 
expected degree of linking. 

3.1 The mean squared writhe of an oriented uniform ran- 
dom walk in a confined space 

In this section, we study the scaling of the writhe of an oriented uniform random 
walk (or polygon) contained in C 3 . 

We are interested in defining the average squared writhe of an n-step uniform 
random walk or polygon in confined space, where the average is taken over the 
entire population of open or closed uniform random walks or polygons in the 
confined space C 3 . We distribute vertices according to the uniform distribution 
on the cube. More explicitly, the space of configurations in this case is O = 
[0, l] 3 ( n+1 ) \ N, and Q = [0, l] 3n \ N respectively, where N is the set of singular 
configurations, i.e. when a walk or polygon intersects itself. Then N is a set of 
measure zero (Randell 1987a, Randell 1987b). 

The average writhe over the space of chains or polygons is zero as there is 
a sign balance occurring due to the mirror reflection involution on the space of 
configurations. This is why we choose to study the mean squared writhe of a 
uniform random walk of n edges in C 3 . 

Theorem 3.2. The mean squared writhe of an oriented uniform random walk, 
or polygon ofn edges, each contained in C 3 , is of the order of 0(n 2 ). Similar 
results hold if C 3 is replaced by a symmetric convex set in M. 3 . 
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Let us consider two (independent) oriented random edges l\ and l 2 of an 
oriented uniform random polygon P n and a fixed projection plane defined by a 
normal vector ( e S 2 . Since the end points of the edges are independent and 
are uniformly distributed in C 3 ,the probability that the projections of l\ and 
I2 intersect each other is a positive number which we will call 2p. We define 
a random variable e in the following way: e — if the projection of l\ and 
I2 have no intersection, e = — 1 if the projection of h and l 2 has a negative 
intersection and e = 1 if the projection of l\ and l 2 has a positive intersection. 
Note that, in the case the projections of l\ and l 2 intersect, e is the sign of their 
crossing. Since the end points of the edges are independent and are uniformly 
distributed in C 3 , we then see that P(e = 1) = P(e = -1) = p, E[e] = and 
Var(e) = E[e 2 } = 2p. 

We will need the following lemma, modeled after the Lemma 1 by Arsuaga 
et al (2007), concerning the case when there are four edges involved (some of 
them may be identical or they may have a common end point): h,l 2 , l[ and 1' 2 . 
Let ei be the random number e defined above between li and l[ and let e 2 be 
the random number defined between l 2 and 1' 2 . 

Lemma 3.3. (1) If the end points ofli,l 2 ,l[ and 1' 2 are distinct, 

then E[eie 2 ] = 0(this is the case when there are eight independent random 
points involved). 

(2) If h = l 2 , and the end points of h,l'i and V 2 are distinct(this reduces 
the case to where there are only three random edges with six independent points 
involved), then E[eie 2 ] = 0. 

(3) If h and l[, or I2 and V 2 have a common end point, then E\e\t 2 \ = 0. 

(4) In the case where l\ — l 2 , the endpoints of l\ and l[ and l\ and 1' 2 are 
distinct, and l[ and 1' 2 share a common point (so there are only five independent 
random points involved in this case), let E[eie 2 ] = u. In the case where l\ and 
l 2 share a common point,the endpoints of h and l[ and l\ and 1' 2 are distinct, 
and l[ and 1' 2 also share a common point ( so there are four edges defined by six 
independent random points involved in this case), let E[e\e 2 ] — v. Finally let 
E[eie 2 ] — w in the case where l\, l 2 , l[ and V 2 are consecutive (so in this case, 
there re four edges defined by five independent random points). Then we have 
q' = 3p + 2(2u + v + w) > 0, where p is as defined before. 

Note that in comparison with Lemma 1 by Arsuaga et al (2007), we have 
included the case in which some of the four or three edges involved are consec- 
utive. 

Proof. (1) This is true since ei and e 2 are independent random variables in this 
case. 

(2) For each configuration in which the projections of l[ and 1' 2 both intersect 
the projection of l\ (since otherwise 6162 = 0), there are eight different ways to 
assign the orientations to the edges. Four of them yield e\e 2 = —1 and four of 
them yield £162 = 1. Since the joint density function of the vertices involved is 
simply , where V is the volume of the confined space C 3 , thus by a symmetry 
argument we have E[e\e 2 ] = 0. 

(3) This is true since in that case e\ = or e 2 = 0. 

(4) Consider the case when the polygon has six edges. Let be the cross- 
ing sign number e between the edges k and lj. Consider the variance of the 
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summation "^2 ^2 (the summation indices are taken modulo 6): 

KK6 i>* 
3 y»-i,<,»+i 



V 



E E 



= E 



l<z<6 j>i 



( 



\ 



E E 

l<i<6 j>i , 

E E + E 

l<z<6 j>i 1<^<6 j>i 

2 E E ^ [ £ y e (j+i)j] 

l<i<6 i>i 

j^i-l,i,i+l,i+2 

2 E E [ £ y' e (i+i)(j+i)] + S [ e '(j+l) £ (»+l)i]) 

l<i<6 J>« 

j^i— 2,»— l,i,i+l,i+2 

2 E [ £ i(i+2) e (i+l)(i+3)] 



(4) 



l<i<6 

Since the e^'s are identical random variables, i.e., they have the same dis- 
tributions, each term in the first summation of the right hand side in the above 
yields 2p, each term in the second summation yields u (see Lemma 1), each 
term in the third and fourth summation yields v and each term in the fifth 
summation yields w. There are 9 terms in the first summation, there are 12 
terms in the second summation, 6 terms in the third summation, 6 terms in the 
fourth summation and 6 terms in the fifth summation. This leads to 



V 



E E 



l<i<6 j>i 



9-2p + 2(l2-u + 6-v + 6-w) 



6{3p + 2(2u + v + w)). 



(5) 



Since V | J2 J2 £i? | > 0, this implies that 3p + 2 (2u + v + w) > 

l<i<6 j>i 

0, as claimed. □ 



Proof of Theorem \3. e A Let us consider the orthogonal projection of P n to a 
plane perpendicular to a random vector £ G S 2 . Throughout the proof the 
averaging, E, is always over the space of configurations. The average, over the 
space of configurations, squared writhe of the projection of P n to that plane is 
given by 
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E [Wrj(P n )] = E 



E E 



= E E +2 E E 



l<i<n j>i 



l<i<n j>i 

j'^2 — 2,2— 1,2,2+1 



+ 2 E E ^N^+i),] 

l<2<n j>£ 

j/i— l,i,i+l,i+2 



+ 2 E E 



(£ + ^ [ e i(j+i) e (i+i)i]) 



l<z<n j>i 

j^i— 2,t— l,i,t+l,i+2 



+ 2 E [ e »(i+2) e (i+l)(i+3)] 
l<i<n 

= n 2 (p + 2(u + v)) - n(3p + 2(4u + 5v - w)) 



(6) 



where p,u,v,w defined as in Lemma 13.31 By Arsuaga et al (2007) it has been 



proved that p + 2(u + v) = q > 0, thus E Wr 2 (P n 
Lemma [33] we can sec that E 



- qn 2 + O (n). Using 

Wr 2 (P n ) is bounded from below by qn 2 — 6qn. 

Let us now take a partition of the surface of the 2-sphere A = {I\ , I2, . . . , I m } 
such that the writhe of the projection of P n is constant in each < j < m. By 
definition, we have that, for a sequence of partitions such that /i (Afc) — > 0, 
the mean squared writhe of P n is equal to 



E [Wr 2 (P n )] = E 



16tt 2 
1 

16tt 2 
1 



E 



E 



(7- lim E Wrt.(P n )8S\ 
lun £ Wr u {P n )8s\ 

l<s<m k I 



lim V Wr 2 E (P n )SS 2 + 2 V Wr CiI (P„) W^ t (P„) <5S 2 

m(a„)-k) \ * i< Sjt < mfc 



lim E 



167T 2 M (A fc )^0 

1 



E ^rf s (P„)<5S 2 + 2 ^ Wr i3 {P n )Wr it {P n )5S 2 



16tt 2 



\ l<s<mfc 1 < s , t < m fc 

(7) 
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where we use the theorem of Lebesgue's dominated convergence, since the func- 
tions Sk (P n ) = (j2i< s <m k Wr^ B (P n ) SS^j , Sk ■ — > K are measurable func- 
tions, bounded above by (^i<s<m ^ r £a (-Pn) ^1 > where Cr^ s (P„) is the num- 
ber of crossings of the projection of P n to the plane perpendicular to £ s and 
(Ei< s < mfc Cr is (P n )5S) 2 <(24n 2 ) 2 . 

The second term in ([7]) is equal to zero, because 



E[Wr s {P n )}=E 



E E 

l<i<n j>i 



= E E S M = °- ( 8 ) 



l<i<n j>i 



But we proved that E 
V£ € S 2 so 



V^r| (P„) = qn 2 + O (n), and E [Wr c (P n )] = 0, 



* - ^ ^ ( E + O (n)) SSA 

\ 1 < s < m. k j 

(gn 2 + O (n)) Qf^ ds) < P [M/r 2 (P„)] < ( 9 n 2 + O (»)) 



dS 



An 



(qn 2 + O (n)) < P (VFr 2 (P„)) < gn 2 + O (n) 



GS 2 



(9) 



Note that in the case of a uniform random walk R n , one has to add to the 
writhe the crossing between the first and last edges, otherwise the proof is the 



same and we obtain the same result, i.e. E 



Wr(R n y Kiqn 2 + 0(n). 



□ 



3.2 The mean squared linking number of two oriented uni- 
form random walks in confined space 

The proof of Theorem 13.21 can be easily adapted in order to provide an analysis 
of the rate of scaling of the mean squared linking number of open chains. Specif- 
ically, we have the following theorem, generalizing Theorem 13.11 by Arsuaga et 
al (2007). 

Theorem 3.4. The mean squared linking number between two oriented uniform 
random walks X and Y of n edges, contained in C 3 , is is of the order of O (n 2 ) . 
Similar results hold if C 3 is replaced by a symmetric convex set in R 3 . 

Proof. For a fixed orthogonal projection of the walks to a plane perpendicular 
to a vector ^ 6 S 2 , adapting Theorem [XT] of Arsuaga et al (2007) to the case of 
open walks, we have E [lk$ (A, Y)] = ^n 2 q + (n) where q > 0. Then, following 
the proof of our Theorem EH we have that E [Lk (A, Y)} = O (n 2 ) . □ 

Note that q = p + 2 (it + v) has the same value in all theorems. By Arsuaga 
et al 2007 it has been estimated to be q — 0.0338 ±0.024. Our numerical results 
confirm this estimation. 
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Figure 2: For all the possible configurations of the edges i, i + 1 and i + 2 such 
that e Vi+2 7^ 0, we have that e M+2 </> i+ i = |<Ai+i| 



3.3 The mean squared self-linking number of an oriented 
uniform random walk or polygon 

The self-linking number was introduced to model two stranded DNA and is 
defined as the linking number between a curve I and a translated image of that 
curve l € at a small distance e, i.e. SI (I) — L (I, l e ). This can be expressed by the 
Gauss integral over [0, 1]* x [0, 1]* = {(#, y) € [0, 1] x [0, 1] \x ^ y} by adding 
to it a correction term, so that it is a topological invariant of closed curves 
(Banchoff 1976), 



Dlj \ l > - 4tt J[0,1]* J[0,1]* I7W-7WI 3 

+ 2ttJ[0,1] | 7 '(t)X7"(*)l 2 ^ ' 

The first term, in the above, is the writhe of the curve which we studied in 
the last section. The second term is the total torsion of the curve, r (I) , divided 
by 2-7T. This measures how much the curve deviates from being planar. The 
torsion of a curve can be expressed as 

r(0= £ Ml) (11) 

l<i<n 

where <j>i{l) is the signed angle between the binormal vectors Bi and Bj+i defined 
by the edges i — + 1 (Banchoff 1976). 

The following theorem concerns the mean squared self-linking number of an 
oriented uniform random walk or polygon. 

Theorem 3.5. The mean squared self linking number of an oriented uniform 
random walk or polygon of n edges, contained in C 3 is of the order 0(n 2 ). 
Similar results hold if C 3 is replaced by a symmetric convex set in M 3 . 

Proof. We will use the definition of the self- linking number given by (fTU|) . i.e. 
SL(l) = Wr(l) + -^ J2i 4>i ■ The proof is based on the fact that the torsion angles 
4>i Vi 7^ 1, n and that the products e,j0i Vj =/= i + 2, k i + 1 are independent 
in the URP model. 

Let P n denote a uniform random polygon in the confined space C 3 . We 
project P n to a fixed plane defined by a normal vector ^ G S 2 . For each pair 
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of edges of the uniform random polygon U and lj we define a random variable 
ejj as we did in the precious section. Then E[e] =0 and E[e 2 } = 2p. For each 
edge we define a random variable fa such that fa is equal to the signed angle 
between Bi and -Bi+i, the normal vectors to the planes defined by the edges 
+ l and i+ l,i + 2 respectively. Then ^ G [— 7r,7r], Vi. Since each vertex of 
the uniform random polygon is chosen with respect to the uniform distribution, 
fa has equal probability of being positive or negative, thus E(fa) = 0, Vi. Now 
let E\4>1] = w and 2?[|0j|] = w'. For each pair of edges fa, 4>j,i ^ 1, j ^ n we 
have that E[fa(f>j] = 0, since fa,4>j are independent random variables in that 
case. We will now compute the mean squared self-linking number 



E 



E ^ + E 1 

l<i<j<n 1 < i < n 

( \ 

E u 

( \ 

E « 

l<i<j<n 



E ^ + 2 E E 

. l<i<n J l<i<j<n l<k<n 



E 



E * 



, Ki<n 



2£; 



E E e «^ fc 

l<i<i<n l<fc<n 



(12) 



In the previous section we proved that 23 
For the second term we have that 



E 



l<i<^<n 



= ]T 2?[0?]+2 ^ Z5[&0 3 -]=tun + 2.E[0i&,]=O(n). 

1 < i < n l^^i^™ 

(13) 

For the third term, we proceed as follows. If j ^ i — 2,i + 2, then eij,(j>k are 
independent random variables for all k, thus 2?[ey0fe] =0. If j = i + 2 then 
£y , </>fc are independent random variables for all k =/= i+l For fc = then there 
are eight different cases that can occur such that E[eu+2fa+i] 7^ 0(see Figure 
[2|). All of them give +|<^_|_i|. Since the vertices of the polygon are chosen with 
respect to the uniform distribution, all the cases have the same probability, thus 
E[eu+2fa+i] = 2?[|</>i+i|] = w' . So finally we have that 

E[Sl 2 {P n )] =qn 2 + 0{n) (14) 

In the case of a uniform random walk R n the self-linking number is not a 
topological invariant and one has to follow a similar averaging procedure as in 
the proof of Theorem 13.21 Finally, the mean squared self-linking number of a 
uniform random walk R n is E [SI 2 (R n )] = 0(n 2 ). □ 



E 



E 
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Remark 3.6. In the case of a uniform random polygon, Theorem 13.51 can 
be proved using the following thinking. For a uniform random polygon P n , 
E [SI 2 (P„)] = E[Lk(P n ,P n J 2 ], where P„ e is the translated image of P n , by a 
distance e, in a direction defined by a random vector. We can then apply the 
same method as for the mean squared linking number of two uniform random 
polygons used in Theorem 13.11 by Arsuaga et al (2007). In Lemma 13.31 by 
Arsuaga et al (2007) we have to add the case of two random edges U in P n and 
the translated edge l\ of P ne such that the one is the translated image of the 
other. But, in that case we have en = 0. This does not change the rate of the 
scaling of the mean squared linking number and we obtain E[Sl 2 (P n )] = 0(n 2 ). 

3.4 The mean absolute value of the linking number of a 
uniform random walk or polygon with a simple closed 
curve in a confined space 

In this section, following the proof of Theorem 4 by Arsuaga et al (2007), we 
analyze the scaling of the absolute value of the linking number between a uniform 
random walk or polygon and a fixed simple closed curve in confined space. 

Theorem 3.7. Let R n (or P n ) denote an oriented uniform random walk (or 
polygon, respectively) of n edges and S a fixed simple closed curve both confined 
in the interior of a symmetric convex set o/R 3 . Then the mean absolute value 
of the linking number between R n (or P n ) and S has a scaling with respect to 
the length of the walk (or the polygon) of the form 



In order to prove this theorem, we will need the following theorem from 
probability theory by C. Stein (Stein 1972). It is used to obtain a bound between 
the distribution of a sum of the terms of an m-dependent sequence of random 
variables (that is X\, X2, . . . ,X 3 is independent of X t , X t +i, . . . provided t~s > 
m) and a standard normal distribution. 

Theorem 3.8. Let x\, X2, ■ ■ ■ , x n be a sequence of stationary and m-dependent 
random variables such that E[xt] — 0, < 00 for each i and 



E[\L{R n ,S)\]^0{^). 



(15) 




) 



2 



< C = lim -E V 

71— >DO 71 \ ' 



< OO, 



(16) 





(17) 



for some constant A > 



12 



Proof of Theorem \3.7\ Note that the confined space can be any convex space 
and a simple closed curve may be of any knot type, but for simplicity we will 
assume that the confined space is the cube given by the set C — {(x,y,z) : 
— \ < X ,V, Z < 5} and that the simple closed curve S is the circle on the xy- 
plane whose equation is x 2 + y 2 = r 2 , where r > is a constant that is less 
than |. Let ej be the sum of the ±l's assigned to the crossings between the 
projections of jth edge lj of P n and S, we need to take the sum since, in this 
case, the projections of lj may have up to two crossings with S. It is easy to see 
that ej = 0, ±1, ±2 for each j, the e/s have the same distributions and that, by 
symmetry, we have E[ej] — for any j. If \i — j\ > 1 mod (n), then e; and £j 
are independent, hence we have E[ei€j] — 0. Let p' — E[e 2 ] and v! — E[ei€i + i]. 
Then, if n = 3, 



= E 



. Ki<3 



Ei 



<i<3 Cl 



= El< l <3^f] + 2 El< lJ <3^] 



3p' + 6u' = 3(p' + 2u') > 



(18) 



Thus, we have p' + 2u' > 0, where p' = E[e 2 ] and u' — E[eiCj]. It follows 



that 



< C = -E 

n 



E 



: p' + 2U' 



(19) 



l<j<n 

for any n. If we ignore the last term e n in the above, then we still have 

2" 



< C = lim 



1 



n — 1 



-E 




= p' + 2u' 



(20) 



Furthermore, the sequence £1,62, ■ ■ ■ , e n -i is a stationary and 2-dependent 
random number sequence since the e^'s have the same distributions, and what 
happens to ei, . . . , ej clearly do not have any affect to what happens to ej+2, ■ ■ ■ , £n-i 
(hence they are independent). 

By Theorem 13.81 there exists a constant A > such that 



P 



= e< - a ] - $ ( a ) 



P[ e ^ ay/(n-l)(j/ + 2u') ] - $(a) 

\l<i<n-l 

P\ £i < w j - $ 



< 



A 



A 

< 



l)(p' + 2u / ) 



< 



(21) 



where w — ay (n — l)(p' + 2m'). 

Note that the linking number between the oriented uniform random polygon 
P n and S is equal to the sum | Ei<j< n e «- 
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Then as n — ► oo, ^ lZi<i<n— l €i — > ^> where Z is a random variable that fol- 
lows the normal distribution with mean and variance a 2 = \{n — l)(p' + 2u'), 



l.C 



N(Q,\(n-l){p' + 2u')) 
lows the half normal distribution and E 

2U')) 1 ' 2 = 0(Vn). 
Thus, 



So the random variable 

2 Sl<i< n -1 e 



4 1 
2 Sl< 



2 Z_jl<i<n— 1 

= i(2/7r(n-l)(p' 



fol- 



E 



KKti-l 



- J5 



2* 



<£[|Lfe(i^,,S)|] < £ 



Ki<n-X 



But J5 



<i<n-l ' 



(22) 



= 0(-y/n) and i£[|e n |] is a constant independent of 

n, so E[\Lk(P n ,S)\] =0(Vn). 

The proof carries through similarly in the case of a uniform random walk 
i?„ and a simple closed curve S. Note that in that case one does not have to 
ignore the last term in (|19|) and has to carry through an averaging procedure 
over all possible projections as well as in the proof of Theorem 13.21 The result 
is again E[\Lk(R n , S)|] = 0{y/n). □ 

Remark 3.9. One can understand the above result, for the case of a uni- 
form random polygon, and a simple closed curve using the following argument 
proposed by De Witt Sumners. We know that Lk(R n , S) equals the algebraic 
number of times the polygon R n passes through the surface Si with S as perime- 
ter. Let to = 0(n) be the number of times the polygon R n passes through the 
surface Si with S as perimeter. We can assume that k — 0(n) of those edges 
are non-consecutive (note that consecutive edges cancel each other and do not 
contribute to the linking number). Then we associate a variable Xi — ±1 to each 
one of those non-consecutive edges depending upon the orientation of the edge. 
Observe that these are independent random variables. Then Xi, 1 < i < k is a 
one-dimensional random walk, thus the distance of the starting point of the ran- 
dom walk of k — 0(n) steps is of the order 0(y/n), thus \Lk(R n , S)| = (9(y / n). 



Remark 3.10. We stress that it is the mean absolute value of the linking num- 
ber of two uniform random walks or polygons in confined space and the mean 
absolute writhe and self-linking number of a uniform random walk or polygon 
in confined space that are of greatest interest for us. Although their analysis 
is much more challenging, we propose that these provide a clearer picture of 
the scaling of the quantities associated with the geometry and topology of the 
chains or polygons. We study these numerically in the next section. It would 
be interesting for future work to prove the numerical scaling analytically. 



4 Numerical Results 



In this section we will describe results obtained by simulations of uniform ran- 
dom walks and polygons in confined space and of equilateral random walks. 
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First we consider the scaling of the mean squared writhe, £[Wr ], and the mean 
absolute value of the writhe, E [|Wr|] of an oriented uniform random walk and 
polygon of n edges in confined space. Then we study the scaling of the mean 
absolute value of the linking number _E[|Lfc|] between an oriented uniform ran- 
dom walk or polygon of n edges and a fixed oriented simple closed curve; and 
the mean absolute value of the linking number between two oriented uniform 
random walks or polygons of n edges each. Finally we study the scaling of 
the mean absolute value of the linking number between two oriented equilateral 
random walks of n edges whose starting points coincide, (ALN), and the scaling 
of the mean absolute value of the self-linking number of an oriented equilateral 
random walk of n edges, with respect to the number of edges, (ASL). 

4.1 Generation of data 

To generate uniform random walks and polygons confined in C 3 , each coordinate 
of a vertex of the uniform random walk was drawn from a uniform distribution 
on [0, 1], and to generate equilateral random walks, each edge vector was drawn 
from a uniform distribution on S 2 . 

For the computation of the linking number or the writhe of uniform or 
equilateral random walks or polygons, we used the algorithm by Klenin and 
Langowski (2000), which is based on the Gauss integral. For each pair of edges 
ei,e2, their linking number is computed as the signed area of two antipodal 
quadrangles defined by the two edges over the area of the 2-sphere. 

We estimated the linking, the writhe and the self-linking numbers between 
oriented uniform random walks and polygons by analyzing pairs of 10 subcollec- 
tions of 500 oriented uniform random walks or polygons ranging from 10 edges to 
100 edges by a step size of 10 edges, for which we calculated the mean and then 
computed the mean of the 10 means for our estimate. We did the same for the 
study of the linking number and the self-linking number of equilateral random 
walks. For the computation of the scaling of the linking number between an 
oriented uniform random walk or polygon and an oriented simple closed planar 
curve, first we considered a fixed square and we analyzed 10 subcollections of 
500 oriented uniform random walks or polygons ranging from 10 edges to 100 
edges by a step size of 10 edges. In order to illustrate that the result holds for 
any fixed knot, we also considered a fixed trefoil and we did the same analysis. 

4.2 Analysis of data 

In this section first we analyze our data on the mean squared writhe of an ori- 
ented uniform random walk and for an oriented uniform random polygon in 
confined space. We fit our data to a function of the form predicted in Theorem 
I3.2l using the method of least-squares. Then we check the variance and the coef- 
ficient of determination. The coefficient of determination takes values between 
and 1 and is an indicator of how well the curve fits the data. We continue 
further our study of the writhe of uniform random walks and polygons by ana- 
lyzing our data on the mean absolute value of the writhe of a uniform random 
walk or polygon in confined space. We do not have an analytical result for this 
scaling, but the data are very well fitted to a linear function. Thus, we observe 
that the distribution of the writhe of a uniform random walk or polygon has 
the property that the scaling of the mean squared writhe is equal to the scaling 
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of the squared mean absolute writhe. We know that if a random variable X 
follows the normal distribution, X ~ N(0, er 2 ), then \X\ follows the half-normal 
distribution and = o^/2pK. Then \X\/a follows the ^-distribution and 

X 2 /a 2 follows the ^-distribution with one degree of freedom and mean 1. Thus, 
E[X 2 /<7 2 ] = 1, E[X 2 ] = a 2 , hence we have that E[X 2 ] « E[\X\] 2 . This indi- 
cates that the writhe of a uniform random walk or polygon in confined space 
follows the normal distribution. 

By Arsuaga et al (2007) it has been proved and confirmed numerically that 
the mean squared linking number of two oriented uniform random polygons in 
confined space has a scaling of the form 0(n 2 ). We stress that in the proof of this 
theorem we have used the fact that the linking number of two closed oriented 
curves is independent of the projection. For two oriented uniform random walks 
(open curves) in confined space we have proved that the mean squared linking 
number has a scaling of the form 0(n 2 ) (Theorem l3.1j) . This is confirmed by our 
data. Since we are interested in the mean absolute value of the linking between 
two uniform random walks or polygons in confined space but do not have an 
analytical result study is challenging, we study this numerically. Our data are 
well fitted to a linear function. This matches the observation of Arsuaga et al 
(2007) that the linking number follows the normal distribution. 

In the same subsection, we study the scaling of the mean absolute linking 
number of a uniform random walk or polygon and a fixed simple closed curve 
in confined space. In Theorem 13.71 we proved that it has a scaling of the form 
0{yfn) and this is confirmed by our data. We note that this is different from 
the mean absolute linking number of two uniform random walks or polygons 
in confined space for which our numerical results indicate a scaling of the form 
0{n). 

Next we analyze numerical data for the mean squared self-linking number 
of uniform random walks and polygons. The data strongly support a scaling of 
the form 0(n 2 ) and, thus, confirm Theorem 13.51 We analyze our data for the 
mean absolute value of the self-linking number of a uniform random walk or 
polygon in confined space. Again, we observe a scaling of the form 0(n) which 
suggests that the self-linking number of a uniform random walk or polygon in 
confined space follows the normal distribution. 

We conclude that all the above measures follow the same type of distribution. 
Furthermore this distribution has the property that the square root of the mean 
squared random variable is equal to the mean absolute random variable. This 
strengthens our intuition that the linking number follows a normal distribution. 

4.2.1 Mean squared writhe and mean absolute writhe of an oriented 
uniform random walk or polygon in a confined space 

Our first numerical study concerns the writhe of a uniform random walk or 
polygon. By Theorem 13.21 the mean squared writhe of an oriented uniform 
random polygon grows at a rate E[Wr 2 ] m qn 2 + 0(n). For comparison with 
this analytical result we calculated the mean squared writhe of an oriented 
uniform random walk of varying length and of an oriented uniform random 
polygon of varying length. 

Results are shown in Figure [3l The black curve in the figure illustrates the 
scaling of the mean squared writhe of a uniform random polygon with respect 
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Figure 3: The mean squared writhe of uniform random walks and polygons. 
Values obtained by computer simulations are shown by triangles and squares 
respectively. The black curve is the graph of the mean squared writhe of a 
uniform random polygon and the dashed curve is the graph of the mean squared 
writhe of a uniform random walk (open chain). 



to its number of edges and is fitted to a function of the form qn 2 + a where q 
is estimated to be 0.0329 ± 0.0002 and a is estimated to be -2.1293 ± 0.9466, 
with a coefficient of determination R 2 = 0.9998. Thus the estimate given in 
the theorem is strongly supported by the data. The dashed curve in the figure 
illustrates the scaling of the mean squared writhe of a uniform random walk 
with respect to its number of edges and is fitted to a function of the form 
qn 2 + a where q is estimated to be 0.0324 ± 0.0002 and a is estimated to be 
-1.1499 ±0.9948, with a coefficient of determination R 2 = 0.9997. Hence, the 
estimate given in the theorem is strongly supported by the data. 

Note that q was estimated to be equal to 0.0329 ±0.0002 and 0.0324 ±0.0002 
respectively, which coincides with the predicted value of q given by Arsuaga et 
al (2007), i.e. q = 0.0338 ± 0.024. 

We next estimate the mean absolute value of the writhe of an oriented uni- 
form random walk or polygon. Results are shown in Figure 0]The black curve 
in the figure represents the mean absolute writhe of a uniform random polygon 
contained in C 3 with respect to its number of edges. The curve is fitted to 
a function of the form a + bn where a is estimated to be —0.2372 ± 0.04262 
and b is estimated to be 0.1460 ± 0.0007, with a coefficient of determination 
R 2 = 0.9998. The dashed curve in the figure represents the mean absolute 
writhe of a uniform random walk contained in C 3 with respect to its number of 
edges. The curve is fitted to a function of the form a + bn where a is estimated to 
be -0.1958 ±0.0630 and b is estimated to be 0.1460 ± 0.0010, with a coefficient 
of determination R 2 = 0.9996. Note that in both cases b w Jq. 
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Figure 4: The mean absolute writhe of uniform random walks and polygons in 
confined space. Values obtained by computer simulations are shown by trian- 
gles and squares respectively. The black curve corresponds to the mean absolute 
value of the writhe of a uniform random polygon and the dashed curve corre- 
sponds to the mean absolute value of the writhe of a uniform random walk (open 
chain) 
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As the number of edges of the polygons increases, we observe a growth at a 
rate 0(n). This suggests the conjecture 

Conjecture 4.1. 

y/E[Wr 2 ] ~ E[VWr 2 } = E[\Wr\] (23) 

Its proof would strengthen our intuition that the mean writhe of a uniform 
random walk or polygon in confined space follows the normal distribution. 

4.2.2 Mean squared self-linking number and mean absolute self- 
linking number of an oriented uniform random walk or polygon 
in a confined space 

Our second numerical study concerns the self-linking of an oriented uniform 
random walk or polygon in confined space. 

The self-linking number is defined as (Banchoff 1976), 

SL(l) = Wr{l) + i-r(l). (24) 

The first term is the writhe of the curve, which we studied in the previous 
subsection and the second term is the total torsion of the curve, r(i) divided by 
2-7T which measures how much the curve deviates from being planar. 

We calculated the mean squared self-linking number of an oriented uniform 
random walk of varying length and of an oriented uniform random polygon of 
varying length. 

Results are shown in Figure [5] The black curve in the figure illustrates 
the mean squared self-linking number of uniform random polygons in confined 
space and is fitted to a function of the form qn 2 + a where q is estimated to be 
0.0331 ± 0.0003 and a is estimated to be -0.0131 ± 1.5226, with a coefficient 
of determination R 2 = 0.9993. The dashed curve in the figure illustrates the 
mean squared self-linking number of uniform random walks in confined space 
and is fitted to a function of the form qn 2 + a where q is estimated to be 
0.0331 ± 0.0002 and a is estimated to be 0.3980 ± 1.1778, with a coefficient 
of determination R 2 — 0.9996. Thus the estimates given in the theorem are 
strongly supported by the data. 

Note that the estimated value of q coincides with the predicted value of q 
given by Arsuaga et al (2007), i.e. q = 0.0338 ± 0.024. 

We next estimate the mean absolute value of the self-linking number of an 
oriented uniform random walk or polygon. Results are shown in Figure [U The 
black curve in the figure illustrates the mean absolute self-linking number of 
uniform random polygons contained in C 3 and is fitted to a curve of the form 
a + bn where a is estimated to be —0.1516 ± 0.1051 and b is estimated to be 
0.1476 ± 0.0017, with a coefficient of determination R 2 = 0.9990. The dashed 
curve in the figure illustrates the mean absolute self-linking number of uniform 
random walks contained in C 3 and is fitted to a curve of the form a + bn where a 
is estimated to be -0.1089±0.0809 and b is estimated to be 0.1463±0.0013, with 
a coefficient of determination R 2 — 0.9994. Note that in both cases b ~ */g.This 
suggests the conjecture: 

Conjecture 4.2. 

y/E[SL 2 ] ~ E[\fSL 2 ] = E[\SL\] (25) 
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Figure 5: The mean squared self- linking number of uniform random walks and 
polygons. Values obtained by computer simulations are shown by triangles and 
squares respectively. The black curve is the graph of the mean squared self- 
linking number of uniform random polygons contained in C 3 . The dashed curve 
is the graph of the mean squared self-linking number of uniform random walks 
contained in C 3 . 
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Figure 6: Mean absolute self-linking number of uniform random walks and 
polygons in confined space. Values obtained by computer simulations are shown 
by triangles and squares respectively. The black curve is the graph of the mean 
absolute self-linking number of uniform random polygons and the dashed curve 
is the graph of the mean absolute self- linking number of uniform random walks. 



Its proof would suggest that the mean self-linking number of a uniform 
random walk or polygon follows the normal distribution. 

4.2.3 Mean squared and mean absolute linking number of oriented 
uniform random walks in a confined space 

Our third numerical study concerns the linking between two oriented uniform 
random walks or polygons in confined space. By Theorem l3.ll the mean squared 
linking number between two oriented uniform random polygons grows at a rate 
0{n 2 ) with regard to the number of edges of the polygons. 

We calculated the mean squared linking number between two oriented uni- 
form random walks of varying lengths and the results are shown in Figure [7] 
The curve in the figure is fitted to a function of the form a + |n 2 where a is 
estimated to be -0.5499 ± 0.9297 and q is estimated to be 0.0346 ±0.0004, with 
a coefficient of determination R 2 = 0.9991 . Clearly, as the number of edges 
of the walks increases, we observe a growth at a rate 0(n 2 ), as expected from 
Theorem I5~l 

Next we calculated the mean absolute linking number between two oriented 
uniform random walks or polygons of varying length and results are shown in 
Figure [U The black curve in the figure illustrates the mean absolute linking 
number of two uniform random polygons contained in C 3 and is fitted to a 
function of the form a + bn where a is estimated to be —0.0300 ±0.0256 and b is 
estimated to be 0.1030±0.0004, with a coefficient of determination R 2 = 0.9999 
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Figure 7: The mean squared linking number of two uniform random walks 
contained in C 3 . 
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Figure 8: The mean absolute linking number of two uniform random walks and 
polygons contained in C 3 . Values obtained by computer simulations are shown 
by triangles and squares respectively. The black curve is the graph of the mean 
absolute linking number of two uniform random polygons and the dashed curve 
is the graph of the mean absolute linking number of two uniform random walks. 
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Figure 9: Mean absolute linking number of a uniform random walk or polygon 
and a fixed square in confined space. Values obtained by computer simulations 
are shown by triangles and squares respectively. The black curve is the graph 
of the mean absolute linking number of a uniform random polygon and a fixed 
square and the dashed curve is the graph of the mean absolute linking number 
of a uniform random walk and a fixed square. 



. The dashed curve in the figure illustrates the mean absolute linking number 
of two uniform random walks contained in C 3 and is fitted to a function of the 
form a + bn where a is estimated to be —0.0251 ± 0.0350 and b is estimated to 
be 0.1044 ± 0.0006, with a coefficient of determination R 2 = 0.9998 . Clearly, 
as the number of edges of the polygons increases, we observe a growth at a rate 
0(n). This suggests the conjecture: 

Conjecture 4.3. 

^E[Lk 2 ] ~ E[VLk*] = E[\Lk\]. (26) 

Remark 4.4. Note that (j2"6f agrees with our intuition and with the numerical 
results obtained by Arsuaga et al (2007) that the linking number of two uniform 
random walks or polygons in confined space follows the normal distribution. 
Indeed, if Lk ~ N(0, a 2 ), where a 2 = 0(n 2 ), then \Lk\ follows the half-normal 
distribution and = a^2/n — 0(n). Then follows the ^-distribution 

and Lk 2 /a 2 follows the % 2 -distribution with one degree of freedom and mean 
1. Thus. E[Lk 2 /<7 2 ] = 1, E[Lk 2 ] = a 2 = 0(n 2 ), hence we have that E[Lk 2 } w 
E[\Lk\] 2 . 

Our next numerical study concerns the linking between an oriented uniform 
random walk or polygon and a fixed simple closed curve in confined space. By 
Theoreml3~7lthc mean absolute value of the linking number between an oriented 
uniform walk or polygon of n edges and a fixed oriented simple closed curve in 
confined space has a scaling of the form 0(y/n). First we consider the oriented 



23 



3.5 h 




3.0 - 



2.5 - 



2.0 - 



1.5 - 



1.0 - 



0.5 - 



n 



20 



40 



60 



NO 



LOO 



Figure 10: The mean absolute linking number of a fixed trefoil and a uniform 
random walk or polygon in confined space. Values obtained by computer sim- 
ulations are shown by triangles and squares respectively. The black curve is 
the graph of the mean absolute linking number of a fixed trefoil and a uniform 
random polygon. The dashed curve is the graph of the mean absolute linking 
number of a fixed trefoil and a uniform random walk. 



square Si defined by the sequence of vertices (0.1, 0.1, 0.5), (0.9, 0.1, 0.5), (0.9, 0.9, 0.5), 
(0.1, 0.9, 0.5), (0.1, 0.1, 0.5) and a uniform random walk or polygon. The results 
of our simulations can be seen in Figure The black curve shows the growth 
rate of the mean absolute linking number of a uniform random polygon and Si . 
The data are fitted to a function of the form a + b^/n, where a is estimated 
to be -0.1665 ± 0.0212 and b is estimated to be 0.1794 ± 0.0029, with a coef- 
ficient of determination R 2 = 0.9980. The dashed curve shows the growth rate 
of the mean absolute linking number of a uniform random walk and S±. The 
data are fitted to a function of the form a + b^/n, where a is estimated to be 
0.0848 ± 0.0153 and b is estimated to be 0.1576 ± 0.0021, with a coefficient of 
determination R 2 = 0.9986. Thus, we can see that the data is consistent with 
Theorem 13.71 

To illustrate that the rate of growth of the mean absolute linking number 
does not depend on the knot type of the fixed simple closed curve, we consider 
the oriented trefoil S2 defined by the sequence of vertices (0.9, 0.5, 0.5), (0.1, 0.5, 0.4), 



(0.5, 0.3, 0.9), (0.6, 0.3, 0.1), (0.2, 0.9, 0.6), (0.5, 0.2, 0.5), (0.9, 0.5, 0.5) and a uni- 



form random walk or polygon. The results of our simulations can be seen in 
Figure 1101 The black curve shows the growth rate of the mean absolute linking 
number of a uniform random polygon and S2. The data are fitted to a function 
of the form a + b^/n, where a is estimated to be —0.1232 ± 0.0381 and b is esti- 
mated to be 0.3610 ± 0.00514, with a coefficient of determination R 2 = 0.9984. 
The dashed curve shows the growth rate of the mean absolute linking number 
of a uniform random walk and S2 ■ The data are fitted to a function of the form 
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a + b^n, where a is estimated to be —0.0142 ± 0.0432 and b is estimated to be 
0.3496 ± 0.0058, with a coefficient of determination R 2 = 0.9978. Hence, the 
data confirm our analytical result in Theorem 13.71 for any fixed simple closed 
curve. 



4.3 Numerical Results on Equilateral Random Walks 

Equilateral random walks are widely used to study the behavior of polymers 
under 0-conditions. It is of great interest to study the scaling of the linking 
number, the self-linking number and the writhe of equilateral random walks and 
polygons. These situations are also much more complex in comparison to the 
uniform models studied in previous sections, since the probability of crossing of 
two edges depends on their distance and the probability of positive or negative 
crossing is independent upon the previous edges. In this section we present 
numerical results concerning equilateral random walks and polygons. It would 
be of great interest to have rigorous proofs of the scaling observed. 

4.3.1 Mean absolute self-linking number of an equilateral random 
walk 



In this section we discuss our numerical results on the scaling of the self-linking 
number of a random walk. 
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The self-linking number of a random walk X is equal to the writhe plus the 
total torsion of the random walk. 

E[\SL{X)\] = E[\Wr(X) + r(X)\] < E[\Wr(X)\] + E[\r(X)\] 

= E[\Wr(X)\] + E[\ £>|] < E[\Wr(X)\] + E?£ |&|] ( 27 ) 

i i 

where E[^2 i \(f>i\] has been proved to be approximately equal to n\ — ^ 
(Plunkett et al 2007, Grosberg 2008). 

Previous numerical results (Orlandini et al 1994) suggest that the average 
of the absolute value of the writhe of an ideal walk increases as y/n, where n is 
the length of the walk. 

Figure fTTI shows the (|^4S'L|) values obtained in numerical simulations of 
ideal random walks. 

We have fitted the computing data points with the function a + b^Jn, leaving 
the two parameters a and b free. Then a was estimated to be —0.3131 ± 0.0313 
and b was estimated to be 0.4014 ± 0.0042, with a coefficient of determination 
R 2 = 0.9991. 



4.3.2 Mean absolute linking number of two equilateral random walks 
whose starting points coincide 

In this section we discuss our numerical results on the scaling of the linking 
number between two equilateral random walks whose starting points coincide. 
Figure [T^] shows the (ALN) values obtained in numerical simulations of ideal 
random walks in a non-constrained linear form. 

By the numerical results presented in the last section, we expect that the 
mean absolute value of the writhe of an oriented equilateral random walk X of 
n steps will have a power-law dependence on the length of the walk: 

(AWR) S3 O(v^). (28) 

Let X = (Xo, X\, . . . , X n ) and Y — (Yo, Y%, . . . , Y n ) denote two oriented 
equilateral random walks of length n whose starting points coincide, i.e. X$ — 
Yq = 0. One can use the scaling of (AWR) to give information concerning the 
scaling of (ALN) as follows : 

We define X — Y to be the oriented equilateral random walk of 2n steps 
(Y n , ...,Yi,Yq = Xo, X\, . . . , X n ). Its writhe then is 

Wr(X -Y) = Wr(X) + Wr(-Y) + 2L(X, -Y) 
^ Wr(X -Y) = Wr(X) + Wr(-Y) - 2L(X, Y) 
=> 2L(X, Y) = -Wr(X —Y) + Wr(X) + Wr{Y) 
=> \L(X,Y)\ < ±(\Wr(X -Y)\ + \Wr{X)\ + \Wr(Y)\) 
=> \L{X,Y)\ < i(0(V2^) + 0(V^) + 0(^)) 

^\L[X,Y)\<0{y/n). (29) 



2G 




We decided therefore to check whether the average of the absolute value of 
the linking number between two ideal walks increases as y/n. We have fitted the 
computing data points with the function a + by/n, leaving the two parameters 
a and b free. Then a was estimated to be 0.0294 ± 0.0130 and b was estimated 
to be 0.0387 ± 0.0018 with a coefficient of determination 0.9839. 



5 Conclusions 



The measurement of the entanglement of open chains is of great interest for 
many applications, such as the study of the properties of polymer melts. In this 
paper, we focused our study in the case of uniform random walks (open chains) 
and polygons in confined volumes. In Theorems 13. 2\ 13.11 13.51 we gave rigorous 
proofs that the scaling of the mean squared linking number, the mean squared 
writhe and the mean squared self-linking number of oriented uniform random 
walks and polygons in confined space, with respect to their length, is of the form 
0(n 2 ). 

Further, we are interested in the mean absolute value of the linking number 
of two uniform random walks or polygons in confined space. In this direction, 
we prove in Theorem 13.71 that the mean absolute value of the linking number of 
an oriented uniform random walk or polygon and a fixed oriented simple closed 
curve in confined space is of the form 0(^/n). Our numerical results confirm 
the analytical prediction, and furthermore suggest that for two oriented uniform 
random walks or polygons in confined space \/E[lk 2 ] ~ 0(n) ~ E[\/lk 2 ]. A 
possible direction for future work would be to prove these results analytically. 
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Ideal random walks are used to model the behaviour of polymers under 6- 
conditions. We have analyzed numerically the scaling of the mean absolute 
value of the linking number between two equilateral random walks of n steps 
and the mean absolute value of the self- linking number of an equilateral random 
walk of n steps. Both appear to scale as 0(\/n). An important direction for 
future work is to complete the analysis and provide proofs for the scaling of the 
self-linking number and that of the linking number of equilateral random walks 
and polygons. 
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